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ABSTRACT 

The physical processes leading to the disappearance of disks around young 
stars are not well understood. A subclass of transitional disks, the so-called cold 
disks with large inner dust holes, provide a crucial laboratory for studying disk 
dissipation processes. IRS 48 has a 30 AU radius hole previously measured from 
dust continuum imaging at 18.7 ^m. Using new optical spectra, we determine 
that IRS 48 is a pre-main sequence AO star. In order to characterize this disk's 
gas distribution, we obtained AO-assisted VLT CRIRES high resolution (R ~ 
100,000) spectra of the CO fundamental rovibrational band at 4.7 micron. All CO 
emission, including that from isotopologues and vibrationally excited molecules, 
is off-source and peaks at 30 AU. The gas is thermally excited to a rotational 
temperature of 260 K and is also strongly UV pumped, showing a vibrational 
excitation temperature of ~5000 K. We model the kinematics and excitation of 
the gas and posit that the CO emission arises from the dust hole wall. Prior 
imaging of UV-excited PAH molecules, usually a gas tracer, within the hole 
makes the large CO hole even more unexpected. 
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Introduction 



Disk around pre-main sequence stars are the birthsites of planets, yet many questions 
remain about the formation processes of planets and the end stages of the disk. Plan- 
ets are expected to leave d i stinctive gaps in the disk as they clear thei r orbital paths 
flArtvmowicz fc Lubowl Il994i Isrvden et all Il999[ ICrida fc Morbidellil boOTf ). Transitional 
disks , in which dust clearing has begun in the inner regions, may be the predict ed objects 



(e.g. iStrom et al.lll989l : ICalvet et al.ll2002l : iBrown et all 120071 iNaiita et all 120071 ). A young 
substellar c andidate recently disc overed in the transitional disk T Cha lends credence to this 
hypothesis (IHuelamo et al.ll201l[ ). However, dust deficits in spectral energy distributions 
(SEDs) can potentially result from other processes including photoevaporation and grain 
growth. Gas observations can help to distinguish between these scenarios. Gas and dust are 
removed near simult aneously in inside-out disk photoevapora tion models and would have 
similar distributions ([Alexander et al.ll2006t lOwen et al.ll201l[ ). At the other extreme, the 
creation of an apparent dust hole through grain growth and settling would cause no deple- 
tion of the gas. Finally, stellar and planetary companions can clear gaps free of both gas 
and dust, but some movement of material across these gaps would be expec ted depending 
on companion size and location ( lArtymowicz fc Lubowlll996l : IZhu et al.ll201l[ ). The creation 
of gaps by companions can en hance other avenues of disk clearing such as photoevaporation 
( lAlexander fc Armitagell2009l ). Gas observations are also of particular relevance in the case 
of young planets as the gas drives both formation and migration. 



CO rovibrational emission is sensitive to warm ' 
from young disks arouri d both Herbig Ae/Be (e.g. 



100-1000 K gas and is c ommonly seen 



Blake fc Booeertl l2004f ) and T Tauri 



stars ( iNajita et al.ll2003l ). These lines are particularly useful as tracer of gas for transition 
disks as the ~1000 K ga s temperatures and high velocities likely imply gas within the dust 
holes ( ISalyk et al.ll2009l ). Spectroastrometric observations of three transition disks confirm 
that the CO fundamental emission generally traces gas within the dust holes but in on e 
case a 7 AU gas hole within a ~30 AU dust hole was found (jPontoppidan et al.l l2008l ) . 
Lack of CO fundamental emission from the inner regions of Herbig Ae/Be stars appears 
slightly n aore common with a handful of disks currently known: HD 141569A (11 AU radius 
gas hole. Icoto et al.lbooel). HP 97084 and HP 100 546 (11 and >8 AU radii respectively, 
van der Plas et al.ll2009l ). iBrittain et al.l (120071 l2009l ) successfully model ed the emission dis- 



tributio ns in both HP 141569A and HP 100546 with UV-fluorescent gas. iPontoppidan et al. 



( l2011al ) finds that it is natural for higher luminosity sources to have larger inner CO-free 
regions, although some sources have even larger inner holes than expected. 



The radiative environment around young high mass stars is harsher than around T Tauri 
stars and allows investigation of disk dissipation processes under different conditions. The 
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majority of transitional disks found to date surround G-type and later stars, although dust 
holes are seen at all stellar masses. Many factors might be responsible for the lower numbers 
including the rarity of local Herbig Ae/Be stars, the potential for a faster transition timescale, 
and difficulty in identifying transition disks around higher mass stars. Herbig Ae/Be SEDs 



have traditionally been classified into almost fiat SEDs 



decline in the far-infrared (group II) (IMeeus et al 



group I) and those with a strong 



2001 



Hillenbrand et al. 



classes represent fiared and self- shadowed disks (IDuUemond fc Doniinikll2004j) . although the 



19921 ). These 



transition between the two is likely continuous (jMeijer et al.l 120081 ). The high luminosities 
of these stars are capable of producing puffed- up inner rims which shadow the disk behind 
( IDuUemond fc Dominikll2004j : lAcke et al.ll2009[ ) and can produce SEDs qualitatively similar 
to transition disk SEDs due to the difficulty in distinguishing between shadowed and miss- 
ing material. However, dust holes are present in disks around at least some higher mass 
sources with millimeter imaging of MWC 758 (A8) an d HD 135344B (F4) clearly resolving 



dust cavities (Andrews et al. 



201 ll : iBrown et al.l 120091 ). Exoplanets have also been directly 



imaged around old er high mass stars at distances that could explain large inner holes (e.g. 



Marois et al.ll2008l ) implying that their younger precursors may be transition disks. 



The Herbig Ae star IRS 48 in p Ophiuchi star formation region has an unusual disk. It 
was first noted for its strong PAH emission and implied strong UV field. Spatially resolved 
mid-IR images of the dust reveal a ring-like structure at 18.7 /xm with a distanc e of 110 AU 



between the brightness peaks and an inner dust hole with a radius of ^3 AU ( IGeers et al. 



20071 ) . A strong east-west directionality is seen leading iGeers et al.l ( 120071 ) to conclude a disk 



inclination of 48° ±8°. Strong polycyclic aromatic hydrocarbon (PAH) molecular emission, 
likely tracing gaseous molecular material, was found to be centered on the star and fill in this 
large grain ring. The CO gas might be expected to have a similar distribution to the PAHs 
and be present inside the dust hole. The VLT-CRIRES spectrograph offers the combination 
of spectral and spatial resolution to resolve the warm gas in this disk. In this paper, we 
present CRIRES observations of IRS 48 showing spatially extended CO rovibration emission 
and examine both the kinematics and excitation to determine the gas structure and origin. 



2. Observations 
2.1. CRIRES 

We observed Oph IRS 48 (WLY 2-48, a = 16'^27'"37M9, S = -24°30'34'.'8, J2000) with 
CRIRES in the M band (4.63-4.87 /im) covering the CO fundamental rovibrational band 



as par t of a large CRIRES program studying a sample of ~ 100 stars (jPontoppidan et al. 



2011bl ). CRIRES is a cryogenic echelle spectrograph behind the Multi- Application Curvature 
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Adaptive Optics (MACAO) system on the 8m VLT 1, Antu flKaeufl et al.ll2004i ). A slit 
width of 0'.'2 was used to provide a resolving power of R ~ 100,000 (A v = 3 km s~^). The 
observations were taken on several different dates in 2007 and 2008 (see Tabled]). 



The data were reduced following the procedures of iPontoppidan et al.l (120081 ) . The data 
were taken with an ABBA 10" dithered nodding pattern permitting a first order correction 
of the infrared background by pair subtraction. Integration times for each position were 1 
minute leading to 4 minute nod cycles. The number of nod cycles in each setting is listed 
in Table [H The CRIRES pixel scale was 86 mas which was subdivided during reduction 
into 43 mas/pixel. The 2-D spectra were linearized and then extracted. The spatially 
extended emission required customized extraction. For an overall spectrum, an aperture of 
1'.'2 was used to include all extended emission. Spectra were also extracted in 2 pixel (86 
mas) apertures to produce a spectral map. All spectra were wavelength calibrated using the 
strong telluric features in the standard star spectra. The standard stars all have spectral 
types of O or B, chosen for their lack of spectral features in the M band. Finally, the many 
weak telluric features were removed by dividing the source spectra by the standard star 
spectra. Remaining strong telluric features with transmission of <20% were blanked from 
the analyzed spectra. 

The majority of the observations were take n with the s li t at a position angle of 90°, 
following the east-west directionality seen in the iGeers et al.l (120071 ) images. Spectra at 0° 
and 180° were also taken to examine the orthogonal distribution of the CO emission (see 
Table [I]). 

Flux calibration was done by normalizing to Spitzer IRAC fluxes. Calibration directly 
from the M-band standard star spectra suffers from variable throughput effects, but in this 
case, direct flux calibration was within 20% of the Spitzer IRAC fluxes. We therefore have 
reasonable confidence that the 5 /im continuum flux has not changed significantly since the 
Spitzer observations. 



2.2. Spectral type of IRS 48 



Luhman & Rieke ( 


1999) 


In contrast. 


Geers et al. 



(120071 ) argued for MO based on an optical spectrum obtained with a fiber-fed multi-object 
spectrograph. However, telluric H2O absorption in the optical spectrum may have been 
mistaken for the TiO band absorption characteristic of MO stars. 



To resolve this discrepancy, we obtained a flux-calibrated, low-resolution {R ~ 1000) 
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optical spectrum using LRIS-ADC on Keck I flOke et al.lll995r ) on 2008 May 28. Although 
the observation spans 3000-9300 A, the star is only detected longward of 5500 A because 
of extinction. Som e weak Ha emission is tentatively detected in our spectrum and in the 
Geers et al.l ( 120071) spectrum, likely indic ating the presence of accretion and consequently gas 
in the inner disk (IMuzerolle et al.ll2004j ). Pfund-/3 emission is also detected in the CRIRES 
spectra, another probable indicator of accretion. The presence and relative strength of 
absorption in the H Paschen and Ca II IR triplet lines clear ly indicate that IRS 48 is an 
A- type star, consiste nt with the B r-7 absorption detected by iLuhman fc Riekd (119991 ). By 
comparison with the iPickled (Il998l ) stellar spectral flux library, we assign a spectral type of 
AOl^ to IRS 48, which corresponds to an effective temperature of 9000±550 K (Figure [1]). 

The extinction is estimated by comparing the cont i nuurn slope in our optical spectrum 
to the continuum slopes of ear l y A-ty pe stars in iPickled (119981 ). The extinction curve used is 
that of IWeingartner fc Draind (120011 ) for Ry = 5.5, which approximates t he opacity of large 



dust grains that are thought to be typical of dense star-forming regions (llndebetouw et al. 
2005! ) . We measure Ay = 11.50±0.25 mag, confirming the presence of substantial intervening 
material in the line of sight, as is common in Ophiuchus. 



Based on the spectral type, extincti on, a J-band magni tude of 10.57 (ISkrutskie et al. 



2006! ) . an d a distanc e of 12 pc to p Oph (ILoinard et al.ll2008l ). the stellar luminosity is 14.3 
Lq. The ISiess et al.l (I2OOOI ) evolutionary tracks then imply that IRS 48 is a ~ 2.0 Mq star 
with an age of ~ 15 Myr at the start of the main sequence. The spectral type and luminosity 
are inconsistent with the 1 Myr age that is typically measured for Ophiuchus. At 1 Myr, a 
star with spectral type AO would have a luminosity of ~ 200 Lq and a mass of ~ 4.5M0. 
Potential sources of error in the calculated luminosity could be a high gray extinction and 
scattering off an edge-on disk, the dust opacity (-Rv'=5.5) used in the extinction correction 
or the distance, which would likely exclude membership in Ophiuchus. We discuss this 
discrepancy further in Section El including additional constraints from the CRIRES data. 



3. Analysis 

Both CO emission and absorption are detected from IRS 48 (Figure [2]). CO emission 
is seen in the v=2-l and very weak v=3-2 lines, while the isotopologues, C^^O and C^^O, 
are in absorption. The ^^CO and ^'^CO v=l-0 lines contain a blend of the two components 
with absorption at the line center and emission at higher velocities. The absorption lines 
disappear at energy levels higher than P(6), indicating that the absorbing gas is cold with 
an excitation temperature of <30 K. The emission lines are all well resolved spectrally with 
CRIRES and are double peaked, including more excited transitions free from absorption. 
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Weak CO ice absorption produces a broad shallow feature between 4.66 and 4.685 and 
the HI Pf/3 line is observed at 4.653 fim, underlying the CO v=2-l R(8) line. The H2 v=0 
S(9) line at 4.6947 /im is not detected. 

The CRIRES AO system produes a small enough PSF in the spatial direction that bright 
extended emission can be seen directly in the sky-subtracted 2-D images. The continuum 
was subtracted based on the mean continuum spatial profile and a scaled telluric standard 
was divided out at every spatial position. The IRS 48 emission lines are spatially extended in 
the 2-D spectra before extraction, particularly clearly in the observations taken at a position 
angle of 90°. Figure [3] shows the ^^CO extended emission in the continuum subtracted 2-D 
spectral trace. 

Extended emission is unusual among the larger sample of mostly T Tauri stars studied 
with CRIRES, whe re most sources show no detectable extended emission in the 2-D images 



(IBrown et al.ll201ll . in prep.). The lines are narrow with a width at line base of 18 km s~^ 
while most of the parent sample show line widths well over 50 km s~^ unless the disk is 
known to be face-on. The narrow width, and thus lack of high velocity gas, is consistent 
with gas motions from large radii and limits the inner extent of the gas, assuming Keplerian 
motion. Figure S] shows the spatial cross section of the clean ^^CO stacked lines. 

By extracting the spectra on a per spatial pixel basis, the physical location of the 
different components is probed (see Figure E])- The CO lines on the central stellar position 
appear in absorption in lines with J < 7 and have little emission in all transitions. Strong 
^^CO v=l-0 and v=2-l emission lines can be seen at distances further from the star. Two 
components are visible: a redshifted component to the west and a blueshifted component to 
the east. The maximum flux is seen at ~30 AU (6 pixels) from the central star. The spectra 
taken at a position angle of 0° show emission in both directions centered around km s~^ 
(Figure El), in agreement with the 90° position angle of the ring seen in the 18.7 /xm dust 



image ( IGeers et al.l 120071 ). There are, however, some asymmetries apparent in the PA=0° 
data with the north side being brighter and more strongly peaked towards the line center. 
A similar brightness asy mmetry between t he north and south sides of the ring is seen in the 



18.7 fim dust (Figure El iGeersetaLl 120071) 



One additional calculation that can be made from the 2-D data is an inclination estimate 
based on the assumption of a symmetrical tilted disk. The major axis is determined to be 
56 AU ± 2 AU peak to peak and 120 AU ± 5 AU in full extent. The minor axis is 42 AU 
± 2 AU peak to peak and 80 AU ± 5 AU in full extent. Assuming a position angle of 90°, 
the inclination is then 42° ± 6°. 
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4. Velocity structure 

4.1. Spectroastrometry 

Spectroastrometry, determining the line flux spatial mean and standard deviation in 
every velocity bin, provides a c onvenient framework to model the emission in IRS 48 (e.g. 



Pontoppidan et al.ll2008l . l2011al ). Collapsing the spatial information in this fashion provides 
several advantages. The biggest source of uncertainty in the direct imaging is almost certainly 
the continuum subtraction. Spectroastrometry, however, is ideally suited for recovery of 
weak contrast line signal from strong continuum. While this is less of a problem in the 
brightest flux spectral bins of the ^^CO lines, the contrast in the line wings and weaker 
transitions of CO v=2-l and ^'^CO is signiflcantly worse. The more limited information also 
makes modeling easier than full 2-D modeling while retaining the main characteristics of the 
spatial distribution. 

The spectroastrometric signal is the flux-weighted center of the spectral trace. The 
weighted mean (x) was found according to the formula x = ^XiFi/UFi where Xi is the spatial 
position and Fj is the flux at that point. The continuum has a dilution effect on the observed 
line signal. The signal that is actually measured has a contribution from both the continuum 
and the line so XTotai{i')FTotai = XTotai{i'){Fc + Fl) = Xc{v)Fc + Xl{i^)Fl, where subscript 
C denotes the continuum component and subscript L represents the line component of both 
the flux (F) and the weighted mean (x). We assume that the continuum flux is pointlike 
and thus xc = 0. The actual line spatial center is therefore xl{^){^ + Fciy) / Fl{i^))i where 
Fc/Fi is the continuum to line flux ratio. The ^^CO lines in IRS 48 are relatively faint 
so Fc/Fl is ~9 in the line center, resulting in signiflcant dilution in the measured signal. 
Fc/Fl becomes increasingly large in the line wings as the signal becomes fainter. 

To measure the width of the emitting region, we turn to the second moment or variance 
(cr^), representing the standard deviation [a) of the signal. It can be found using the formula 
cr^ = Tj{xi—PYFi/TjFi where Xi and Fi are as above. P is the central spatial position taken as 
the peak position of the continuum emission at 4.7 /im. The variance is similarly continuum 
diluted such that (j'^otaiW)FTotai = (^TotaMiFc + Fl) = crl{u)Fc + (tI{p)Fl. However, 
the continuum does have a standard deviation due to the telescope PSF and seeing so 
al ^ 0. Here, al = a'^otaiO- + Fc{v)/Fl{v)) - al;Fc{v) / Fl{v). A further correction to 
examine the line width around the line center, so P = xl-, can be made by subtracting x| 
so a\ p^^^ = '^L ^ xi- The plots present which can then be compared with the PSF 
determined from the width of the continuum signal. Directly calculating the variance in the 
line center gives an emitting region width of 4 AU t}^^ centered around the line center. 
However, this estimate is likely too large as each spectral pixel covers a range of radii due to 
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the curvature of the ring. A more detailed and accurate discussion of the hne extent using 
simple models can be found in Section 14.21 

Figure E] shows the detected fluxes, spatial centers and variances of the ^^CO, ^^CO 
v=2-l and ^^CO emission lines. The spectroastrometric signals are calculated over 25 pixels 
(129 AU) which includes all the flux, with larger ranges only increasing the noise. Lines are 
stacked to enhance the signal. Only lines without line overlap and central absorption are 
included. For ^^CO, this includes P(7), P(8) and P(ll). For ^^CO, it includes R(9), R(ll), 
R(12) and R(13). CO v=2-l has the most clean lines as even the low energy rotational levels 
are free of absorption, including P(l), P(2), P(3), R(l), R(2), R(5),and R(6). The line shapes 
of all the isotopologues are similar but the absolute measured signals are different due to 
stronger continuum dilution in the weaker lines. Multiplying by the correction factors (2.8 for 
CO v=2-l/C0 v=l-0 and 4.6 for ^^CO/CO v=l-0), results in identical spectroastrometric 
signatures within the noise (Figure |6]). All of the transitions are therefore emitting from the 
same physical region. 

Figure [7| shows the spectroastrometric signatures arising from the different epochs. Ob- 
servation in three different epochs 3-9 months apart allows us to search for variability. While 
the continuum width varied between the different nights due to AO performance, the CO 
structure appears stable over the year timescale. 



4.2. Kinematic modeling 

A simple Keple rian disk model is fitted to the data to examine the physical structure 



of the CO (following iPontoppidan et al.ll2008l ). The model is gridded in the disk frame with 



rings of gas at the Keplerian velocity. Each ring is assigned a temperature according to an 
overall disk temperature profile. The surface density is assumed to be constant with a total 
column density, Nqo- The local line width is set to 3 km s~^, the instrumental spectral 
resolution. The fiux at each grid point is determined for the different rotational lines. The 
grid is then transformed into the observer's Cartesian coordinates by twisting to the correct 
position angle and inclination. The emission is so extended that the 0'.'2 slit does not capture 
all the light, so the data is corrected for the truncation accounting for the AO-corrected PSF 
of 0'.'18 based on the measured continuum width (see Fig. [7]). The emission in the PA=0° 
position is strongly limited by this parameter. 

A grid of models varying the stellar mass, inclination, inner radius {Rq) and position 
angle were run to examine the interdependence of these parameters. The two most interesting 
parameters are the inner radius and the stellar mass. In order to robustly determine the 
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inner radius regardless of mass, the minimum value for each mass was determined. Figure 
|9] plots mass against Rq, assuming the inclination associated with the minimum value. The 
hole size (Rq) minimum is 30 AU throughout the mass range with optimal inclinations of 
50° at the low mass end and 30° at the high mass end (Figure [H]). 

Stellar mass is degenerate with inclination in these models. Accurately determining the 
incli nation is the r efore the best method of determining the stellar mass. The inclination 



from lGeers et al.l (120071 ) is 48° ± 8°, which becomes incompatible with masses higher than 3 
Mq. The CRIRES 2D images give a slightly more face-on inclination measurement of 42°±6° 
(see Figure [TO] shows the values for stellar mass versus inclination with the CRIRES 
inclination overlaid. The stellar mass of 2 Mq ± 0.5 Mq agrees with the determination from 
the optical spectra and is incompatible with a 1 Myr AO star. 

As the surface density is held constant, the temperature profile is also a proxy for 
the distribution of gas. We test two different distributions, a power-law and a uniform 
temperature ring. The power law distribution is of the form, T = To{R/Ro)^'^, where i?o 
is taken to be the inner most edge of the gas and the power-law declines from there until 
the outer edge of the disk, set at 100 AU. Tq is set at 260 K in agreement with the rotation 
diagram results from section ISTTl A steep power law index (g>0.4 for best results) is required 
in order to show the double peaked flux profile. The uniform temperature ring is centered 
at Ro and Tq but has width AR. The effects of varying emission widths in a uniform 
temperature ring with constant density are shown in Figure [Hi The width of the emitting 
region is directly related to the variance, al- The best fits result from the narrowest emitting 
regions, limited only by the grid cell size. Overall there was little difference between the two 
distributions as the good fits from both concentrated the emission close to Rq with a width 
of a few AU or less. Physically this narrow region is likely the dust hole cavity wall. 

In summary, the best fit to the data was obtained using the derived stellar mass of 2 
Mq (g2]) and an inclination of 42° from the 2-D CRIRES data (Figure [8]). The hole size is 
30 AU with the disk at a position angle of 95°, in good agreement with the mid-IR image. 
The model fit of the PA=0° observations is not as good due to the asymmetries in the data 
(i.e. Figure E]). The structure appears more complex than a simple uniform ring in that 
direction. 
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5. Excitation 

5.1. Rotational temperatures 

Rotation diagrams, which use the relative strengths of the rotation lines to understand 
the gas temperature and column density, were constructed for both the emission and absorp- 
tion components. Interpretation of rotation diagrams is commonly based on the Boltzmann 
distribution and the assumption of LTE level populations, so 

m = !^^{Ej/kT)^ 

gj Quot 

where n is the number, Quot is the partition function, and Ej is the energy of level J. We 
plot In (Fj/gAu) against i^upper/^ where Fj is the line flux, g is the degeneracy (2J[/p+l), A 
is the Einstein A coefficient, 9 is wavenumber in cm~^ and -Eupper is the upper state energy. 
In the optically thin case, a linear fit will match the data well. The slope (a) is related to 
the temperature, Troj such that Tjiot = — 1/a. The intercept, y, is related to the CO column 
density, Nco, such that 

Nco = Qnote'^^ (2) 
where d is the distance and Aemit is emitting area. 



5.1.1. Absorption lines 



The absorption lines are strongest in the lowest energy rotational lines and absorption 
features can be seen in the less abundant isotopes, ^^CO, C^^O and C^^O. The absorption 
lines are spectrally unresolved at 3 km s~^ and have a central vt,sr of 3.8 km s~ ^ in agreement 
with the p Ophiuchi region (e.g. iDi Francesco et al.ll2004t iMaruta et al.ll2010l ). The derived 
gas excitation temperatures from the different isotopologues are ~25 K (see bottom half of 
Table |3]). The ^^CO absorption is also highly optically thick. These cold temperatures are 
consistent with an origin in foreground cloud material. The absorption lines are therefore 
probably not related to the disk and are excluded from further analysis. 



5.1.2. Emission lines 

Quantifying the emission line fluxes is complicated by the contaminating absorption 
lines at low J. A clean template of the emission line proflle is made by stacking the higher 
J lines which are free from absorption. This template proves to be a good match for the 
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CO v=2-l and ^^CO emission lines as well as the CO v=l-0. The profile is then fitted 
via chi-squared minimization on top of a measured local continuum to derive each line flux. 
This method provides more consistent measurements than integrating the flux directly which 
is strongly affected by noise. The low J ^^CO lines are contaminated at low velocities by 
absorption, so a fit was made to the line wings. 

The rotation diagram (Figure [TT!) for ^^CO v=l-0, ^^CO v=2-l and ^^CO relative to the 
energy of J=0 in each vibrational level indicates that the different transitions have similar 
properties. The distributions all appear linear, in agreement with optically thin gas at a 
single temperature. In this case, the derived rotation temperatures are ~250 K (Table [3]). 
Small changes in the derived slope (i.e. temperature) can result in large changes in derived 
population. All the slopes are in the range 100-260 K, but the better constrained species 
are at the higher end. The total mass {Nco^^emit) is well constrained, but the column 
density is not due to the lack of constraints on the emitting area. Based on the kinematics 
( §4.2p . we estimate that the area is a ring at 30 AU with a width of 0.5 AU and use this area 
to calculate the column densities in Table |3l A gas temperature of 250 K is slightly high 
compared to the hkely dust temperature of ~150 K at 30 AU based on the low luminosity of 
the star. However, the gas temperature may be much higher than the dust temperature due 
to strong UV radiation. The CO v=2-l lines are too strong to be purely thermally excited 
with a vibrational temperature much higher than the rotational temperature. We attribute 
the excess excitation to U V fluorescent excitation which is discussed in Section 15.21 The low 
CO v=l-0 lines may be slightly optically thick but the errors are large due to the central 
absorption lines and a linear flt is well within the errors. However, the optically thin flts 
produce a low isotopic ratio ^^CO/^^CO of 15, challenging this assumption. 

In order to determine a limit on the amount of CO within the hole, flux limits were 
determined from the 2-D images. In general, the flux levels were similar to residuals from the 
continuum subtraction in regions without CO emission and were on the order of 10% of the 
peak flux. Assuming optically thin gas, rotation diagrams for gas at different temperatures 
were compared with the flux limits to determine the maximum CO gas mass that could 
be within the hole without detection (Figure [T2|) . The amount is strongly dependent on 
temperature with much lower limits at high temperatures and weaker constraints at low 
temperatures. In general, the inner region should be hotter than the gas further out. The 
mass limit is ~10^° g for temperatures above 500 K. This translates to a gas density of ~10^ 
cm~^, assuming a uniform distribution within the 30 AU hole and a disk height resulting 
from hydrostatic equilibrium. 
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5.2. UV fluorescent excitation 

The distance of the gas from the star, resulting in cooler temperatures, and the pres- 
ence of high energy CO v=2-l and v=3-2 lines suggest that a non-equilibrium process such 
as UV fluorescence is likely responsible for their line emission. The measured vibrational 
temperature is 2200 K, greatly in excess of the thermal rotational temperature of 260 K. A 
strong UV field is to be expected from an A type star such as IRS 48 and the strong PAH 
emission from within the hole requires its presence. 

A simple UV fluorescence model was run to examine the likelihood of the stellar UV 
field producin g the observed vibrationally excited emission at 30 AU. The model follows 



Brittain et al.l ( 120071 ) but expands the vibrational levels covered to 35 X^S^ and 25 A^U 
levels. The model assumes statistical equilibrium so each level maintains a steady popula- 
tion. We then considered excitation from X^S"*" to A^U (qa-x) balanced by spontaneous 
emission (Aa-x) back to the ground electronic state. Decay out of the excited electronic 
state populates a range of vibrational levels in the ground state, including vibrational levels 
which are too energetic to be thermally populated. We also include spontaneous relaxation 
between vibrational levels in the ground electronic state (Ax-x), which occurs primarily 
between levels with Av = 1. The transitions covered between the ground state X^E"*" elec- 
tronic state and the A^U excited electronic state occur at ultraviolet wavelengths 1270 to 
4000 A, but the most populous levels are mainly driven in the range 1270 to 1700 A. The 
model also includes coUisional excitation, assumed to be thermalized at Tq = 260 K. 

In general, the model output depends purely on the input UV spectrum. While accre- 
tion shocks produce much of the UV emission in lower mass T Tauri stars, the similarity 
between the photospheric and shock temperatur es in Herbig Ae stars results in the photo- 



spheric emission dominating (jValenti et al.ll2000l ). Based on our derived stellar properties 



we constructed a 9000 K blackbody and compared this to an HST STIS UV spectrum of AB 
Aur (from StarCAll^. AB Aur is a good proxy for IRS 48 with a spectral type of AO and 
age of 4 Myr. While the spectrum contains many lines, the overall flux and shape is a good 
match to the 9000 K blackbody, once scaled to the appropriate luminosity. The calculations 
used the observed spectrum in the available range (1140-1729 A), which accounts for 99% of 
the transitions in this source, and the blackbody at longer wavelengths. 

The observed fractional population in each vibrational level was calculated using the 
rotational diagrams (Figure [TTl) . The energies of each rotational state were normalized to 
the zero point of each vibrational level. The intercept was used to derive the population in 



^http://casa. Colorado. edu/^ ayres/StarCAT/ 
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each level, which was then divided by the total population. The ground state v=0 is by far 
the most populated. The low signal-to-noise in many of the faint lines introduces a lot of 
scatter in the rotational diagram. Given the similarity in line profiles and spectroastrometry, 
the gas is likely all close in temperature so we fix the rotational temperatures at 260 K and 
then derive the populations. 

Figure [13] shows the vibrational diagram with the overlaid model fit. A linear fit to 
the data (v > 0) gives a vibration temperature of 2200 K, significantly higher than the 
rotational temperature. However, the thermal component to the v'=l level is significant, 
actually causing the vibrational temperature to be underestimated. The model of the effect 
of the UV field at 30 AU, including the thermal contribution from gas at 260 K under 
the assumption of LTE excitation, is shown in red while the purely vibrational excitation 
is shown in blue. For v=2 and 3, the blue points lie completely under the red as there is 
virtually no collisional excitation to these levels. The vibrational temperature of the model is 
~5000 K. The model fit is in excellent agreement with the data, indicating that the expected 
UV field from IRS 48 for = 2 M0 and L^, = 14.3 Lq can produce the observed energy 
distribution. 

We also investigate the expected excitation from the younger, higher mass star implied 
by the evolutionary age of p Oph ( §2.2p . The luminosity is set at 200 L© and the distance 
is kept at the 120 pc distance to p Oph. The black crosses in Figure [13] represent the 
expected excitation from such a star at 30 AU. While the vibrationally excited emission 
increases by about an order of magnitude, the collisional excitation also increases due to 
a higher temperature of 500 K at 30 AU. It appears unlikely that IRS 48 is a younger, 
more luminous, but highly extincted star. One remaining potential source of error in the 
luminosity would be an incorrect distance, implying that IRS 48 is not actually in the p Oph 
molecular cloud. The vibrational excitation does not provide a constraint on this scenario 
as a larger distance simply corresponds to a large hole size and the vibrational diagram is 
unchanged. 



6. Discussion 



Based on our CRIRES data, the CO emission around IRS 48 comes from a ring at 30 
AU in Keplerian rotation around the central star. We propose that the CO gas is thermally 
and UV excited in th e dust wall a t the h ole edge. The radius of the gas emission and the dust 
hole edge seen in the lGeers et al.l (120071 ) 18.7 pm dust image are similar, indicating that the 
CO and dust rim are co-spatial. Our Keplerian model is consistent with a ring of emission 
with little to no flux required from outside the hole wall region, as seen in the required steep 
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power law index. The majority of the stellar UV light must reach the hole wall in order to 
excite the CO v=2-l and v=3-2 emission seen along with the v=l-0 emission, given the good 
match between the fluorescence model and data. While little dust shielding is likely present 
inside the hole, the PAH molecules could absorb some of the UV field. Although without 
knowing either the abundance or distribution of the PAHs, it is difficult to estimate the 
magnitude of the absorption. The hole wall may be vertically puffed up so that a fraction 
of the stellar UV light reaches the rim unattenuated. 

Despite the lack of CO interior to 30 AU, so me material i s prese nt in the inner disk. 
The strong, centrally-peaked PAH emission seen by iGeers et al.l (120071 ) indicates that small 
particles remain within the dust hole in IRS 48. Generally, PAHs are small enough that they 
are expected to be dynamically coupled to the gas. The presence of PAHs also promotes 



the formation o: 



(IJonkheid et al. 



H2, CO and other small molecules even in the absence of small grains 



2OO6I). O t her si gns of material inside the hole are the Ha lines in both our 



2008 and the iGeers et al.l (120071 ) optical spectra and the Pf/3 line present in all the CRIRES 
data - all probable signs of ongoing accretion. We know little about the distribution of gas 
within the hole from these signatures although material is likely still reaching the star. 

Any overall picture of the IRS 48 system must therefore explain both the truncation of 
the dust grains and CO at 30 AU as well as the destruction of CO, but not PAHs and atomic 
hydrogen, within the hole. The strong UV field within the hole may be the cause of the 
different survival rates. Larger PAH molecules have a bette r chan ce of surviving in high UV 
environments than smaller ones. According to IVisser et al.l ( 120071 ). PAHs of sizes 50 carbon 
atoms and smaller are destroyed at the r=l surface out to 100 AU around a Herbig Ae/Be 
star. Given the likely low opacity of the disk interior only PAHs with 100 carbon atoms or 
greater would be able to populate the hole outside of about 5 AU, sufficient to explain the 
mid-IR VISIR images. On the other hand, CO has a photodissociation lifetime of < 1 year 
inside the hole without shielding. Preliminary modeling indicates that gas densities below 
10^ cm~^ are inadequate for rapid reformation of CO and prevent self-shielding (Bruderer, 
priv. comm.). Our limits on the CO gas density within the hole lie well below this value 
for temperatures greater than 500 K. However, the CRIRES data can give little constraint 
on any low temperature gas. Gas with temperatures of 300 K or less could be present at 
densities of 10^ cm~^. This regime of PAH emission but no CO emission places stringent 
limits on the amount of gas present in the hole. 

The sharp truncation of CO and dust points towards either physical truncation by a 
companion or photoevaporation of the CO. While the UV field around IRS 48 is strong and 
the accretion signatures are weak, it is not clear that p hotoevaporative inside-out clearing 



would leave the material inside the hole in this state (e.g. [Alexander et al.ll2006l lOwen et al. 



- 15 - 



20111 ). Gas and dust inside the hole in the photoevaporative clearing scenario is usually 
residual material still draining onto the star and is not being replenished. This system 
would thus have to be in a specific time period where the gap has open, the density of the 
interior material has dropped below the CO photodissociation threshhold but has not yet 
drained completely. A companion, on the other hand, would naturally reduce the fiow of 
material inwards and create a drop in density of both gas and dust. A final possibility, 
although probably the least likely, is that no material is fiowing through the hole and the 
PAHs are formed in situ via the collision of large rocky bodies. 

The CO hole in the Oph IRS 48 disk is the largest gas ring discovered to date for any 
pre- main sequence star. The hole is more than twice as large as the 11 AU holes around 
the Herbig Ae/Be stars, HD 141569 and HD 97084, which have similar spectral types to 
IRS 48. SR 21, a G-type star, with a ring of emission at 7 AU is the lowest mass source 



with a known large near-IR CO hole (jPontoppidan et al.l 120081 ). Vibrational emission of 



similar extent to the v=l-0 lines in all cases indicates that the UV field may be the driving 
force in maintaining the presence of the CO holes via photodissociation. Large near-IR CO 
holes of >1 AU have not yet been found in K- or M-type T Tauri transition disks, despite 
the presence of substantial dust holes. In fact, many transitio n disks contain CO close 



to the star at rad ii significantly smaller than the dust hole size (jPontoppidan et al.l 12008 



Salyk et al.l l201l[ ). The UV energy from T Tauri stars is significantly smaller than from 



higher mass stars with a large portion generated in accretion shocks. If transiti on disks have 



even smaller accretion rates than classical T Tauri stars as has been suggested ( iNajita et al. 



20071 ). T Tauri stars are unlikely to have both an absence of absorbing inner dust and a 
strong UV field, hampering the creation of large photodissociated CO holes. Although CO 
fundamental emission is detected from the inner dust wall of IRS 48, any non-detection of 
CO emission from T Tauri disks with holes (e.g. CoKu Tau/4) may be caused by difficulty 
exciting the 4.7 fim CO lines at large radii due to both the weaker UV field and lower 
temperatures, rather than a lack of gas at the inner rim of the dust disk. Further work to 
better understand the excitation of gas in holes is needed to determine if the creation of such 
large CO holes is a unique phenomenon of intermediate mass stars. 

We are also left with an evolutionary puzzle regarding IRS 48. The CRIRES data, both 
from kinematics and excitation, are much more consistent with a lower mass star of around 
2 Mq and a luminosity of ~15 Lq, as derived from the optical spectroscopy. The data are 
inconsistent with a prototypical 1 Myr AO star with a mass of 4.5 Mq and luminosity of 200 
Lq. However, the resulting age from evolutionary t racks of ~15 Myr is much older than the 



1 Myr age of Ophiuchus (ILuhman fc Riekelll999l ). with the star being underluminous by 



more than an order of magnitude. Such an old PMS star might be expected to have a debris 
disk of reprocessed dust and no gas, but this is clearly not the case for IRS 48. IRS 48 appears 
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to be a member of Ophiuchus with proper motions in agreement with the cloud (E. Mamajek, 
priv. comm.), a position close to the core, emission lines centered at t he cloud veloc it y and 
an extinction consistent with a location within or behind the cloud. iBaraffe et al.l (|2009[ ) 
suggest that episodic cold accretion leads to different contraction timescales, and thus some 
stars may appear much older from evol utionary mode l s than their true age. While IRS 48 is 
more massive than the stars studied by lBaraffe et al.l (120091 ). similar mechanisms may apply 
even at higher mass and help explain the apparent age discrepancy for IRS 48. If the stellar 
parameters, cloud membership and PMS evolutionary tracks are all correct, Oph IRS 48 
would be an extremely old member of Ophiuchus, possibly left over from a previous epoch 
of star formation. 



7. Conclusions 



We have used the combination of high spatial and spectral resolution of CRIRES to 
observe spatially extended CO emission from Oph IRS 48. IRS 48 was known to be a 
transitional disk from the ima ged 18 Atm dust r ing but the centrally peaked PAHs lead to 
an expectation of interior gas (iGeers et al.ll2007| ). Despite the presence of PAHs within the 
dust hole, the CO emission peaks at 30 AU, co-spatial with the dust hole edge, and shows 
no emission interior to the dust hole. The spatial distribution of the vibrationally excited 
gas and isotopologu es is the sanie, in dicating a common origin. The CO geometry is in 
agreement with the iGeers et al.l (120071 ) dust image with an east-west major axis, brighter 
emission from the north side compared to the south, and a derived inclination of 42°±6°. 
The gas is thermally excited to a temperature of 260 K, but is also strongly fluoresced by 
the central star's UV field, showing vibrational excitation temperatures of ~5000 K. The 
new spectral type of AOl° resolves th e discrepancy of strong PAHs around what was thought 
to be an MO star (IGeers et al.l 120071 ) and explains the vibrationally excited emission. Our 
kinematic modeling is consistent, both spectrally and spatially, with a ring of gas at 30 AU 
in Keplerian rotation with a width of only a few AU or less. All the evidence points towards 
the gas being UV excited in the outer rim of the hole. A sudden drop in density to <10^ 
cm~^, as might be caused by a companion, could cause the CO to photodissociate on short 
timescales but maintain the larger PAH population. IRS 48 is therefore an excellent target 
for planet finding follow-up. 
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Fig. 1. — Comparison of the H Paschen series in the IRS 48 LRIS-ADC spectrum (black) 
with the iPickled (119981 ) spectral library. The presence and strength of these lines clearly 
indicates that IRS 48 is an A-type star, with a best fit of AOl^. 
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Fig. 2. — VLT-CRIRES spectrum of IRS 48 with the upper row showing close ups of different 
hnes and the lower panels showing the emission lines in detail on the top and the full depth of 
the absorption lines on the bottom. ^^CO/'^CO and C^^O are seen in absorption, particularly 
strong in low energy J transitions. The ^^CO v=l-0 and v=2-l lines show strong double 
peaked emission features. 
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Fig. 3. — Continuum subtracted 2D spectra showing the extended hne emission in the ^^CO 
P(7) and P(8) hues from 2007 September 5. The emission peaks around 30 AU and is much 
weaker towards the central position. 




Fig. 4. — Cross section of the stacked lines at PA=90°. The black line is the total spatial 
distribution while the blue and the red lines are the blue and red sides of the line. The 
peak emission averaged over the entire line width occurs at ± ~25 AU. This is slightly 
smaller than the largest peak to peak extents, marked with dashed vertical lines, due to 
the inclusion of gas not at apoapsis of the inclined ring. When compared directly to the 
continuum distribution, the CO spatial profiles are consistent with a small emitting region 
(see also §4.2). 
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Fig. 5. — (Left) Spatial changes in the profile of the P(8) transition of ^^CO v=l-0 in 
increments of 86 mas (10 AU) for PA=90° (left) and PA=0° (right). Emission occurs from 
the outer regions while line of sight absorption occurs towards the central star at lower J 
(<6). Keplerian rotation can clearly be seen in the emission with the emission peak in the 
top spectra redshifted compared to the bottom. The missing data on the left, removed due 
to contamination by the telluric P(8) line, show the clear separatio n between science and sky 
lines. (Right) The slit positions overlaid on the iGeers et al.l (120071 1 18.7 /im image. Some of 
the asymmetries of the dust are echoed in the gas, including brighter emission in the eastern 
peak compared to the west and brighter emission from the north compared to the south. 
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Fig. 6. — Comparison of the flux (top), spatial center (middle) and extent (bottom) of the 
IRS 48 CO hnes for ^^CO (left), CO v=2-l (middle), and ^^CO (right) at a position angle of 
90°. The v=l-0 lines have been scaled to the appropriate flux ratio and overp lotted as dotted 
red lines in the v=2-l and ^^CO panels. The bottom panel shows the standard deviation 
of the line emission. The standard deviation of the continuum signal is shown as a dotted 
horizontal line in the ^^CO panel. Deconvolving the PSF results in an emitting region width 
of 4 AU ^4Au^- The agreement in spatial postion and line shape of the different transitions 
indicates that this emission comes from the same physical region. 
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Fig. 7. — Comparison of the flux (top), spatial center {xTotah middle top), and seeing mea- 
sured from the continuum {(TTotaU bottom) of the IRS 48 CO lines from different observation 
dates. Black and blue are 2 Aug, 2008 at PA of 90 and 270 respectively (spatial center for 
270° is multiplied by -1). Red is 5 September, 2007 and green is 3 May 2008. 
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Fig. 8. — Comparison of the flux (top), spatial center (middle) and extent (bottom) of the 
IRS 48 CO lines. The left column is the spectroastrometry at position angles of 90° and 270° 
while the right column is 0° and 180°. The overplotted red line is a simple model of a ring 
of gas in Keplerian rotation at 30 AU around a 2 Mq star with an inclination of 42°. The 
temperature of the disk outside this radius is taken to have a power-law distribution with 
an index of -0.6. 



- 29 - 




1.5 2.0 2.5 3.0 3.5 4.0 4.5 60 80 100 120 

Stellar moss (Mg) Position Angle 



Fig. 9. — (left) plot of Rq at each M, choosing the inclination and position angle which 
produce the best fits. The inclinations range from 30° to 50° with the higher masses tending 
to better fits at smaller inclinations. Rq is therefore fairly robust across a range of stellar 
masses, (right) The effects of varying emission region widths in a uniform temperature ring 
with constant density. The emission is clearly dominated by a small region close to i?o- 
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Fig. 10. — plot of stellar mass vs inclination. The white star marks where the measured 
inclination of 42° is at the minimum. The red vertical error bars show the error on the 
inclination while the diagonal line traces the lowest values within the inclination errors. 
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Fig. 11. — Rotation diagram of the CO emission lines. Filled circles are ^ CO, diamonds 
are ^^CO (dotted fit), stars are ^^CO v=2-l, and triangles are ^^CO v=3-2 (dashed fit). The 
energies refer to those of the rotational levels relative to J=0 in the respective vibrational 
state. 
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Fig. 12. — Upper limits on CO mass within the 30 AU hole assuming different gas temper- 
atures. 
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Fig. 13. — Vibrational diagram for IRS 48. The black points are the observed fractional 
population in each vibrational level. The linear fit indicates a vibrational temperature of 
2200 K, greatly in excess of the rotation temperature of 260 K. The red circles are the 
expected vibrational fractions from a UV fluorescence model of IRS 48, including the thermal 
population. The blue circles, underlying the red for v>l, are the vibrational component in 
each level before the collisionally excited population is added. The black crosses represent 
the expected excitation from our 1 Myr old, high mass spectral parameters with a luminosity 
of 200 Lq at 30 AU. While the vibrationally excited emission increases by about an order 
of magnitude, the bulk of the increase is from the expected increase in coUisional excitation 
with the increased temperature at 30 AU. 
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Table 1. Summary of CRIRES observations 



Date 


Wavelength 


Position 


Number of 


Standard 




(nm) 


Angle 


Nod Cycles 


Star 


2007 Sep 5 


4730 


90 


1 


HR 5812 


2008 May 3 


4730, 4833 


90 


2 


HR 6084 


2008 Aug 2 


4710 





3 


HR 6084 


2008 Aug 2 


4730 


0,90,180,270 


3 


HR 6084 


2008 Aug 5 


4868 


90 


3.5 


HR 6175 



Table 2. Model parameters 



Parameter 


Best fit model 


Range 


Hole size (Rq) 


30 AU 


25-35 AU 


Stellar mass 


2 Mo 


1.5-2.5 M0 


Inclination 


42° 


36-48° 


Position angle 


95° 


85°-100° 


q 


0.6 


>0.4 


AR 


1 AU 


< 3 AU 
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Table 3. Summary of rotation diagrams^ 



Molecule 


Rotational 


CO Mass 


CO Column Density 




Temperature (K) 


(g) 


(cm-2) 


12C0 v=l-0 


256 


7.2 10^3 


7.2 10^^ 


i^CO v=l-0 


267 


4.6 10^2 


4.6 10^^ 


i^CO v=2-l 


204 






i^CO v=l-0 


26 






i^CO v=l-0 


17 






Ci^O v=l-0 


30 







^Emission rotation diagram parameters are above the horizontal 
divide and absorption parameters are below. 



